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Abstract 

The dependence of the pseudoscalar meson masses and decay constants on sea 
and valence quark masses is compared to next-to-leading order (NLO) Chiral Per- 
turbation Theory (ChPT). The numerical simulations with two light dynamical 
quark flavors are performed with the Wilson-quark lattice action at gauge coupling 
/3 = 5.1 and hopping parameters k = 0.176, 0.1765, 0.177 on a 16^ lattice. 0{a) 
lattice artifacts are taken into account by applying chiral perturbation theory for the 
Wilson lattice action. The values of the relevant combinations of Gasser-Leutwyler 
constants L4, L5, Lq and Lg are estimated. 

1 Introduction 

The low energy dynamics of strong interactions in the pseudo-Goldstone boson sec- 
tor of QCD is constrained by the non-linear realization of spontaneously broken 
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chiral symmetry [T]. In an expansion in powers of momenta and light quark masses 
a few low energy constants - the Gasser-Leutwyler constants - appear which pa- 
rameterize the strength of interactions in the low energy chiral Lagrangian [2^ . The 
Gasser-Leutwyler constants are free parameters which can be constrained by ana- 
lyzing experimental data. In the framework of lattice regularization they can be 
determined from first principles by numerical simulations. In experiments one can 
investigate processes with different momenta but the quark masses are, of course, 
fixed by Nature. In numerical simulations there is, in principle, much more freedom 
because, besides the possibility of changing momenta, one can also freely change 
the masses of the quarks. This allows for a precise determination of the Gasser- 
Leutwyler constants - once the simulations reach high precision. First steps towards 
this goal have recently been done by several authors [21 1^ El El including our Col- 
laboration iziiHiini. 

The main difficulty for numerical simulations in lattice QCD is to reach the 
regime of light quark masses where ChPT is applicable. The reason is the critical 
slowing down of simulation algorithms for small quark masses and lattice spacings. 
We apply the two-step multi-boson (TSMB) algorithm ^01 which allows to perform 
simulations with small quark masses within the range of applicability of next-to- 
leading order (NLO) ChPT |71IH1. 

Another important aspect of investigating the quark mass dependence in numer- 
ical simulations is the possibility to use ChPT for the extrapolation of the results to 
the physical values of u- and d-quark masses which would be very difficult to reach 
otherwise. In fact, ChPT can be extended by changing the valence quark masses in 
quark propagators independently from the sea quark masses in virtual quark loops 
which are represented in the path integral by the quark determinant. In this way one 
arrives at Partially Quenched Chiral Perturbation Theory (PQChPT) (11111211131 ■ 
The freedom of changing valence and sea quark masses substantially contributes 
to the power of lattice QCD both in performing quark mass extrapolations and in 
determining the values of the Gasser-Leutwyler constants |14j . 

For a fast convergence of numerical results to the continuum limit it is important 
to explicitly deal with the leading 0{a) lattice artifacts. An often used method is 
the application of the 0{a) improved lattice action ^Hl- We apply an alternative 
technique which in the pseudo-Goldstone boson sector is equivalent to the 0{a) 
improvement of the lattice action. In this method the (unimproved) Wilson action 
is used in the Monte Carlo generation of gauge configurations and the 0{a) effects 
are compensated in PQChPT itself. This means that we apply chiral perturbation 
theory for the Wilson lattice action. Our calculations showed that in practice this 
method gives results with good precision |^. 
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The plan of this paper is as follows: in the next two sections we collect the NLO 
(PQ)ChPT formulas for ratios of pseudoscalar meson masses and decay constants. 
In section 21 a discussion of the general form of the NNLO tree-graph corrections is 
also included. In section 0] the results of numerical simulations is presented. The 
last section is devoted to a summary and discussion. 



2 Valence quark mass dependence 

In this paper we use the notations introduced in [0] which slightly differ from those 
of ref. Ul] and JB]- The dimensionless variables for quark masses and 0{a) lattice 
artifacts are denoted, respectively, by 

2Bomg 2Woacsw 

XA = —72— ' PA = 72 ■ (1) 

Jo JO 

Here mg is the quark mass, a the lattice spacing, Bq and Wq are parameters of 
dimension mass and (mass)^, respectively, which appear in the leading order (LO) 
chiral effective Lagrangian, csw is the coefficient of the 0{a) chiral symmetry break- 
ing term and /o is the value of the pion decay constant at zero quark mass. (Its 
normalization is such that the physical value is /o — 93MeV.) For fixed sea quark 
mass xs the dependence of the pseudoscalar meson mass and decay constant on the 
valence quark mass xv can be described by the variables 

c - XV _ PS /„x 

C = — , VS = — • (2) 

XS XS 

For instance, in case of a number of Ng equal mass sea quarks the ratios of decay 
constants are given by 

Rfvv = ^ = 1 + 4(e - l)xsLs5 

JSS 

~ 64^^ + ? + 2??5) log + ^^(1 + Vs) log(l + Vs) , (3) 

and similarly 

Rfys^ffs^^^ 2(e - i)xsLs, + 7T^(e - 1) - 77^(1 + log 

JSS GANsTT'^ GANsTT^ 1 + r/5 

"128^^ + C + 2r]s) log + -^^(1 + vs) log(l + r]s) . (4) 

I^Sk {k = 5) denotes the relevant Gasser-Leutwyler coefficient at the scale fo^/Xs■ 
This is related to Lk defined at the scale /o and defined at the generic scale /j, 
according to 

Lsk = Lk-Ck logixs) = L'k - Ck log(-|xs) • (5) 



with the constants Cfc, {k = 4,5,6,8) given below. Similarly, the corresponding 
relations for the coefficients Wsk introduced in jl6j are: 

Wsk = Wk-dk log(xs) = Wl-dk log(4xs) • (6) 

The constants in ^ and @ are given by 

1 Ns _ (jVK2) _ (iVf-4) 

^ 256^2 ' C5 256^2 ' C6 5i2A^|7r2 ' 512A^,7r2 ' ^' ' 

respectively, 

^ " 2567r2 ' " 2567r2 ' " 256iV|7r2 ' " 256A^,7r2 ' ^""^ 

For the valence quark mass dependence of the (squared) pseudoscalar meson 
masses one can consider, similarly to (jSl and (jlj), the ratios 

Rmyy ^ , Rmvs = ^ . (9) 

In the present paper we prefer to divide these ratios by the tree level dependences 
and consider 



Rnvv = = 1 ~ 7 — 



and 



16A^,7r2 ^ 



2m?.,, (£ - 1) 

+4(C - l)xs{2Ls8 - Ls5) + 8Ns%^vsXs{Ls4 - Wse) 

(4 + ij 

^(l + 2??5)log(l + r/s) 



XS (g^ + g + ??5 + 3T?gg) 
16iV,7r2 + 1) 



log(e + r]s) . (11) 
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A useful quantity is the double ratio of decay constants .17, which does not 
depend on any of the NLO coefficients. In other words there one can see the chiral 
logarithms alone. The NLO expansion for this quantity is: 



fvs 1 , XS' /t i\ XS , ^1 ^ + VS 
f f = 1 + ooAr 2 ^ - 1 - ooAr 2 ^ + log — 

JVVJSS 32iV^7r^ 32A's7r^ 1 + Vs 



The double ratio of the pion mass squares |18j corresponding to (fTUIl and Ullf) 
has the NLO expansion 

''''^ - (e + i?miymis e(e + 1) 

XsiC"^ + ^ + VS+ 3??5g^) log(e + ris) _ xs{2r)s + l)log{l+Vs) 
16A^sVr2^(C + 1) WN^tt"^ 

im:^^ — + — a^TT) - ^''^ ■ 

2.1 Quadratic corrections 

A complete NNLO (i. e. two-loop) calculation in PQChPT for our physical quantities 
is not yet available. Nevertheless, the general form of NNLO tree-graph ("counter- 
term insertion") contributions can be given JH]^- For instance, one has for the pion 
mass square: 

StTI ad 2 / \ / \9 /2 2\ / A\ 

2 = <^IXS + 0^2X3 [XA + XB) +a3{XA + XB) + OiiiXA + Xb) ■ (14) 

^AB 

(Here A and B denote generic quark indices: S is the label for the sea quarks, V 
for valence quarks.) For the pion decay constant there is a similar expression. This 
information is very useful in order to estimate the importance of the NNLO terms 
in our present range of quark masses. 

The general characteristics of the NNLO terms is that they are proportional to 
the quark mass square: x's- (Here we only consider terms in the continuum limit 
and hence neglect lattice artifacts. This will be to some extent justified a posteriori 
by the observed smallness of 0{a) terms.) Neglecting loop contributions, which are 
at the NLO order relatively small, the dependence on the quark mass ratio ^ is at 
most quadratic and can, therefore, be represented by terms proportional to (^ — 1) 
and (^ — 1)^. Therefore, these contributions have the generic form 

DxXsiC-l) + QxXsi^-l? ■ (15) 



^We thank the authors for communicating us the content of this paper prior to publication. 
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Here X denotes an index specifying the considered ratio as, for instance, X = 
fVV, nVS etc. for the single ratios and X = fd and X = nd for the double 
ratios RRf and RRn, respectively. The NLO tree-graph contributions for the single 
ratios Rf and Rn are also proportional to — 1). These can be parametrized as 
LxXs{£, — 1) (foi' instance, we have Lfyy = 4L5 and Lnvv = 8(2L8 — L^)). The 
inclusion of Dx-type terms is equivalent to a linear dependence of the effective Lx's 
for fixed xs'- 

L^f = Lx+ DxXS . (16) 

At this point one has to remember that mathematically speaking - in order to 
completely remove the effect of higher order terms - Lx is defined in the limit 

The NNLO coefficients are not all independent but satisfy the relations 
DfVS = i^DfVV , Dnvs = -j^DnVV , 

Dfd = , Dnd = ^, 
Qfd = '^QfVS — QfVV + '^^}vV ' Qnd = '^QnVS — QnVV + -^L^VV ■ i^"^) 

The first line is a consequence of the general structure of the NNLO tree-graph 
contributions. The last two lines follow from the definition of RRf and RRn if one 
only considers NLO and NNLO tree-graph contributions. 

We shall see in section |1] that in our range of quark masses the NNLO tree-graph 
contributions of the form (|15j) are important but can be approximately determined 
by global fits. In this way the NLO constants are better determined. Observe 
that a determination of the -Dx's is only possible in our analysis if different sea 
quark masses are included (see below). 

2.2 0{a?) corrections 

The idea of including leading lattice artifacts in the low energy effective Lagrangian 
for the Wilson lattice action can be extended to higher orders in lattice spacing. In- 
deed, in writing this paper we have seen two recent publications about the inclusion 
of 0{a?) corrections |2U1 121j . The general formulas derived in these papers for the 
0(0^) terms imply that in the formulas for the pion mass-squared ratios ()1U() . (jllj) 
and there are only very little changes. In fact, the changes can be summarized 
by the replacement 

Vs{Lsi-Wsi,) ils{Lsi-WsQ)+'^{NMsA + Ws^-2N,Ws&-2W's^) . (18) 
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Here Wg^ and Wgg denote some new low energy constants appearing in the O(a^) 
part of the effective Lagrangian. This means that fitting the valence quark mass 
dependence with our formulas H1U|) . and H13|) effectively takes into account also 
0(0^) corrections. 

Concerning the ratios of the pion decay constants in Q , © and (|12|) the situ- 
ation is expected to be similar but there, in addition to the O(a^) terms, also new 
types of 0{amg) terms may appear. 



3 Sea quark mass dependence 

The dependence on the sea quark mass can be treated similarly to the valence quark 
mass dependence considered in section |2I Here one chooses a "reference value" of 
the sea quark mass XR and determines the ratios of the coupling and decay constant 
as a function of 

a = — , r = — . (19) 

XR PR 



Instead of r one can also use 



which satisfy 



_ PS _ PR i^f.-. 

Vs = — , VR = — (20) 

XS XR 



^ = ^ . (21) 



With this we have for the decay constants 
Rfss = ^ = l + A{a- l)xR{NsLm + Lm) + 4(r?5fT - r]i^)xR{NsWm + Wr^) 

JRR 

-^CT(1 + r?s) log[a(l + Tjs)] + ^(1 + m) iog(i + m) (22) 

and for the mass squares 
2 

Rnss = = l + T]s-m + 8(<T - l)xRi2NsLRe + 2Lr^ - N.Lr^ - Lr^) 

+8(r/sa - m)XR{'iNsWR(i + 2VF^j8 - N^Wr^ - Wr^ - N,Lr^ - Lr^) 

+ TrTW<^ + 2^5) log[a(l + r^s)] - TTTTwi^ + 2r/ij) log(l + Vr) ■ (23) 

Of course, the coefficients Lr^ and WRk {k = 4, 5, 6, 8) are now defined at the scale 
fo\/XR therefore in the relations © and © Xs is replaced by XR- 

The logarithmic dependence of LskS and WskS have to be taken into account 
also in simultaneous fits of the valence quark mass dependence at several sea quark 
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mass values. Choosing a fixed reference sea quark mass XR we have from Q and 
© with 11 = foy/xR 

Lsk = Lru - Ck log a , Wsk = WRk - 4 log a . (24) 

The NLO PQChPT formulas for the valence quark mass dependence in terms of 
the reference sea quark mass are obtained by the following subsitutions in ©, @, 

(Unii-jini): 

XS crxR , Lsk Ljik , Wsk WRk , 

log(l + r/s) log[cr(l + r]s)] , 

log(C + r]s) log[cr(C + r]s)] , 

log(l+C + 2??5) ^ log[(7(l + e + 2r?s)] . (25) 

An important feature of both the valence- and sea-quark mass dependences 
considered in the present work is that they are ratios taken at a fixed value of the 
gauge coupling (/?). These are renormalization group invariants independent from 
the Z-factors of multiplicative renormalization since the Z-factors only depend on 
the gauge coupling and not on the quark mass. Taking ratios of pion mass squares 
and pion decay constants at varying quark masses has, in general, the advantage 
that quark mass independent corrections - for instance of 0{a) and/or ©(a^) - 
cancel. 



4 Numerical simulations 

We performed Monte Carlo simulations with Ns = 2 degenerate sea quarks on 
a 16^ lattice at /3 = 5.1 and three values of k: kq = 0.176, ki = 0.1765 and 
K2 = 0.177. For the reference sea quark mass we choose kr = kq = 0.176. A 
summary of the simulation points is reported in tabled where also the set-up of the 
TSMB algorithm for the different simulation points can be found. The gauge field 
configurations collected for the evaluation of the physical quantities are separated 
by 10 TSMB update cycles consisting out of boson field and gauge field updates 
and noisy correction steps. It turned out that these configurations were statistically 
independent from the point of view of almost all secondary quantities considered. 
Exceptions are ro/a and Mr (see below) where autocorrelation lengths of 2-5 units 
in the configuration sequences appear. 

We investigated for each simulation point the valence quark mass dependence 
of the pseudo-Goldostone boson spectrum and decay constants; the values of the 
valence k considered for each simulation point are reported in table El In these 
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Table 1: Parameters of the simulations: all simulations were done 
at (3 = 5.10 with determinant breakup Nf = 1 + 1. The other 
TSMB-parameters are: the interval of polynomial approximations 
[e, A] and the polynomial orders ?t-i,2,3 (10) . 



run 




configurations 


e 


A 


ni 


^2 


^3 





0.1760 


1811 


4.50 • 10-4 


3.0 


40 


210 


220 


1 


0.1765 


746 


2.50 ■ 10-4 


3.0 


40-44 


280 


260 - 340 


2 


0.1770 


1031 


3.75 • 10-^ 


3.0 


54 


690 


840 



intervals the valence quark masses are approximately changing in the range i^rrigea < 

nivalence ^ ^TTl^ga' 

A rough estimate of the sea quark mass range can be obtained by considering 
the quantity Mr = (rom-j^)'^, which for the strange quark gives ~ 3.1. (Here ro ~ 
0.5 fm is the Sommer scale parameter which characterizes the distance scale intrinsic 
to the gauge field.) In our simulation points the value of ranges between Mr ~ 
2.10 and Mr ~ 1.09, corresponding to about | and | of the value for the strange 
quark mass. Since the valence quark masses roughly go down to m^aience — ^nisea, 
they reach m^alence — \nis. In our configuration samples we did not encounter 
problems with "exceptional gauge configurations" - in spite of the smallness of the 
valence quark mass. This means that the quark determinant effectively suppresses 
such configurations. 

Standard methods for the extraction of the relevant physical quantities have 
been applied (a more detailed description is given in our previous paper [7j and 
in [121 )• Statistical errors have been obtained by the linearization method |231 124j 
which we found more reliable than jack-knifing on bin averages. 

Within a mass-independent renormalization scheme - defined at zero quark mass 
- the Z-factors of multiplicative renormalization depend only on the gauge coupling 

and not on the quark mass (k). Similarly, the lattice spacing a is also a function 
of the gauge coupling alone Therefore, since our simulation points are at fixed 
gauge coupling /? = 5.1, the ratios of the sea quark masses can be obtained by 
taking ratios of the measured bare quark masses in lattice units Zqamq. Here Zq 
is the multiplicative renormalization factor for the quark mass which is the ratio 
of the Z-factors of the pseudoscalar density and axialvector current [Zq = Zp/Za) 
because we determine the quark mass by the PCAC-relation: = mq'-'^'-' 0. 
(Of course, in the valence quark ratios the factor Zqa also cancels trivially.) The 
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Table 2: Values of the valence quark hopping parameter. 



run 





1 


2 




0.1760 


0.1765 


0.1770 


'I'valence 





1685 





1710 


0.1743 







1705 





1718 


0.1747 







1720 





1726 


0.1751 







1730 





1734 


0.1754 







1735 





1742 


0.1759 







1745 





1750 


0.1763 







1750 





1758 


0.1767 







1770 





1772 


0.1775 







1775 





1778 


0.1779 







1785 





1785 


0.1783 







1790 





1791 


0.1787 







1800 





1797 


0.1791 



obtained values of the sea quark mass ratios fjj = mgi/mgo {i = 1,2) are given in 
table 121 together with some other basic quantities. 

Note that by identifying the quark mass ratios in the ChPT formulas with the 
ratios of the PCAC quark masses ( "axialvector Ward identity quark masses") one 
assumes that these two kinds of renormalized quark masses are proportional to each 
other. As it is shown, for instance, by eq. (48) in j21j this is indeed the case - apart 
from lattice artifacts of 0{amg) and 0{a'^). The quark mass independent part of the 
0(0^) terms are cancelled by taking ratios. The remaining quark mass dependent 
lattice artifacts are neglected in the present paper. 

The critical value of the hopping parameter where the quark mass vanishes can 
be estimated by a quadratic extrapolation using the values of ai^2- 

^ _ niqi _ {K' ^ — Kg/) + dcr{K' ^ — K^/)^ ^^g-j 
niqO (Kq Kcr ) ~l~ da{^^Q ^cr )^ 

The values of ai^2 in table |21 give the solution: Kcr = 0.1773(1) and d^j = —11.2(8). 
(The relatively large absolute value of d^ shows that the quadratic term in the 
extrapolation is important.) 

The value of the lattice spacing a can be inferred from the value of r^/a at 
This can also be determined by a quadratic extrapolation of the values of 
ro/a given in tableElwith the result: ro(Kcr)/a = 2.65(7). Taking, by definition, 
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''o('^cr) = 0.5 fm this gives for the lattice spacing: a = 0.189(5) fm. 

The physical volume following from the lattice spacing is comfortably large: 
L ~ 3.0 fm. Since the minimal value of the pion mass in lattice units in our points 
is am^™ ~ 0.43 for sea quarks and am^™ ~ 0.30 for the lightest valence quark, we 
have Lm-jj > 4.8. 

Another information given by the values of Mr is an estimate of the quark mass 
parameter xs in the ChPT formulas. For instance, in the reference point we have 
from ro/o ^ 0.23 (SE]: Xr**™"*^ ~ Mr/{rofof ~ 39.8. 



Table 3: The values of some basic quantities in our simulation 
points. Statistical errors in last digits are given in parentheses. 









/«2 


ro/a 


2.149(15) 


2.171(88) 


2.395(52) 




0.6747(14) 


0.6211(22) 


0.4354(68) 


Mr = {rom^f 


2.103(26) 


1.824(41) 


1.088(47) 


Zqarriq 


0.07472(32) 


0.06247(51) 


0.03087(36) 


o-i = rriqi/mqQ 


1.0 


0.8361(52) 


0.4132(34) 



4.1 Valence quark mass dependence 

For a fixed value of the sea quark mass xs the valence quark mass dependence of 
the ratios Rfwys-, RT^Vvys-, RRf and RRn is determined by five parameters: 

Xs , XsLs5 , XsLs85 = XsC^Lss - Ls5) , XsLs4W6 = Xs{LsA - Ws&) , rjs . 

(27) 

The dependence is linear in the first four of them but it is non- linear in 775. 

After performing such fits of the data we realized that the sea quark mass de- 
pendence is not consistent with the NLO PQChPT formulas. In particular, the best 
fit values of the xs's have ratios considerably closer to 1 than ai^2 in tablelHland the 
change of the L^'s with xs is also not consistent with ()24() . This shows that NNLO 
effects are important and, therefore, we tried fits including NNLO tree-graph terms 
of the form given in (|15|) . The list of the relevant NNLO parameters is: 

XjiDfVV,nVV , XjiQ fVVJV SJd,nVV,nV SM ■ (28) 

Qfd and Qnd have to satisfy the quadratic relations given in the last line of (|17|) 
but in order to keep linearity we did not impose these relations and fitted the eight 
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parameters in (|28j) independently. After performing the fits one can check how well 
the relations for Q/d Qnd ^-re fulfilled. 

The global fit of the valence quark mass dependence for several values of the sea 
quark mass has twelve linear parameters: the first four in (|27|) with xs replaced by 
XR 

XR , XrLr5 , XrLr85 = XrC^L^is- Lr^) , XrLr4W6 = XRiLR4,-Wm) (29) 

and the eight in (|28j) . In addition there are the non-linear parameters, in our case 
three of them: r]s = ??o,i,2- 



Table 4: Values of best fit parameters for the valence quark mass 
dependence. Quantities directly used in the fitting procedure are in 
bold face. 



Xr 


33.5(2.4) 






XrLr4W6 


5.24(38) ■ 10-2 


Lrawg 


1.564(71) ■ 10-3 


XnQnd 


6.5(1.8) ■ 10-3 


Qnd 


5.80(79) ■ 10-6 


XrLr5 


10.06(44) ■ 10-2 




3.00(19) ■ 10-3 


x\Dfvv 


-9.3(1.7) ■ 10-2 


Dfvv 


-8.3(1.9) ■ 10-5 


XnQfvv 


-2.80(19) ■ 10-2 


QfVV 


-2.50(50) ■ 10-5 


XnQfvs 


-2.197(45) ■ 10-2 


QfVS 


-1.96(29) • 10-5 


xlQfd 


-0.99(14) ■ 10-2 


Qfd 


-0.89(45) ■ 10-5 


XrLrs^ 


-2.10(12) ■ 10-2 


Lrs5 


-6.25(52) ■ 10-^ 


X%,DnVV 


-1.67(20) ■ 10-1 


DnVV 


-1.49(10) • 10-^ 


XjlQnVV 


-8.44(67) ■ 10-2 


QnVV 


-7.53(48) ■ 10-5 


XuQnVS 


-4.05(25) ■ 10-2 


QnVS 


-3.61(29) • 10-5 



Multi-parameter linear fits are easy and, except for degenerate situations, the 
chi-square always has a unique well-defined minimum. Non-linear fits involving the 
ry's are more problematic, therefore we adopted the following procedure: performing 
non-linear fits at individual sea quark mass values we obtained the starting values 
of '70,1,2 • Then for fixed values of ?7o,i,2 we performed a linear fit of the twelve 
parameters in H28() - (|29|) and looked for a minimum of the chi-square as a function 
of %,i,2- For the sea quark masses we imposed the relation xs = ^^XR for the 
NLO parameters the relations in (|24|) with the values of (Ti^2 given in table IHI (The 
possible dependence of the NNLO parameters D and Q on a has been neglected.) 
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The minimum of the chi-square after the non-hnear minimization is near 

r)o = 0.07 , rji = 0.03 , r?2 = 0.02 . (30) 

The minimum as a function of ?7o,i,2 is rather shallow but definitely within the 
bounds < ??o,i,2 ^ 0.10. The minimization of the chi-square of the linear fit does 
not change the r/'s substantially: already the starting values are close to ()30() . This 
confirms the small value of r]s found in our previous paper at /? = 4.68 [S]. 

In contrast to the stable values of the r/'s there are large fluctuations in the basic 
parameter XR- one can obtain values in the range 13 < XR ^ 40 depending on 
the set of functions fitted, on the fit interval etc. This is presumably the effect of 
our small number (only three) of sea quark masses. In order to obtain more stable 
results we fixed ?7o,i,2 according to (|30|) and first determined in a linear fit the three 
parameters XRi XrI^R^W^ and xjiQnd from RRn. These parameters were then used 
as an input in the linear fit of the remaining nine parameters. 

All 18 valence quark mass dependences considered can be reasonably well fitted. 
The best fit is shown by figures ^ and l2j The sum of the chi-squares of the linear 
fits is x^ — 300 for a number of degrees of freedom n.d.f. = 18-12 — 12 = 204. Most 
of the chi-squares comes from the points with largest and smallest valence quark 
masses where there are obviously some systematic deviations, too. The parameters 
of best fit are given in table 0] The values in the table show that there are some 
discrepancies in both relations in the last line of ()17() but the deviations are not 
very large. The first and second relation give: -0.89(49) • 10-5 ~ 2.05(39) • 10"^ 
and 0.52(9) • 10"^ ~ 0.92(8) • 10-^ respectively. 



Table 5: Values of combinations ofak's obtained from the best fit values in table^and\B 



as 


2.24(20) 










ass = (2a8 - as) 


0.762(49) 






(a4 - c^e) 


2.36(9) 


a4s = (2a4 - as) 


2.40(26) 


A4//0 


22.9(1.5) 


1^45 


-1.7(1.8) 


«6845 = (4a6 + 2a8 - 2a4 - as) 


0.658(86) 


A3//0 


6.51(57) 


1^6845 


-5.43(60) 



The values of the NLO and NNLO parameters themselves are also shown in the 
right hand part of table 01 with errors determined (as always) by the linearization 
method. With the help of the formulas in ©-© one can also transfer these results 
to the corresponding L's and W's at some other renormalization scale different from 
fo^/XR- Going to the conventional renormalization scale // = 47r/o and multiplying 
by an overall factor 1287r^ one obtains the values of and uj^ shown in table |SJ 
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Due to the unexpected smallness of the 0{a) contributions it is interesting to 
try a hnear fit of the valence quark mass dependences setting all 0{a) terms to 
zero: ijq = iji = r]2 = 0. This is a fit with eleven parameters because in the formulas 
Ls4,w6 is always multiplied by rjs. The result is a reasonable fit but the chi-square 
is by about 10% larger then in the case of t/o,i,2 7^ 0. The best fit values of the main 
parameters are in this case: XR = 36.1(1.0), 05 = 2.08(14), as5 = 0.502(48). 

The NNLO tree-graph contributions are rather important especially at k = 
0.176. From the point of view of the NLO formulas the situation becomes better 
at K = 0.177 but NNLO is still not negligible there: see figure |H1 (At k = 0.1765 
we have, of course, an intermediate situation between k = 0.176 and k = 0.177.) 
In general, the NNLO contributions are more important in the ratios RnVV a-nd 
RnVS than in Rfvv and Rfvs- In fact, the ratios RnW and RnVS at k = 0.176 are 
dominated by NNLO. The relative importance of NNLO terms is stronger for ^ > 1 
than for ^ < 1. In the double ratios RRn and RRf the NNLO terms are relatively 
unimportant. 

4.2 Sea quark mass dependence 

The results from the fit of the valence quark mass dependence can also be used 
in the investigation of the sea quark mass dependence according to ((^ - (P^ . In 
particular, the values (and errors) of XR and ?/o,i,2 are relevant there. Besides these 
values and the known ratios of the sea quark masses ai^2 (see table two extra 
parameter pairs appear, namely, for Ng = 2: 

Lm^ = 2Lm + Lr^ , Wm^ = 2Wm + Wm (31) 

in ((211) and 

LmsA5 = iLn6 + 2LR8-2Lm-LR5 , VFr6845 = iWm+2WRs-2WR4-WR5 (32) 
in 

Since we only have three sea quark mass values and therefore two independent 
values of Rfss and Rnss a "fit" actually means solving for the four unknowns. 
The results are collected in table IHl The corresponding values of the a's and w's 
are contained in table IHl In this table also the values of the universal low energy 
scales A3^4 are given. (For the definitions see |271 126j or eq. (10) in ^.) Once the 
parameters -L/J45 and L_r6845 are known it is possible to extrapolate the continuum 
NLO curves (without the 0{a) contributions) for Rfss and Rnss to zero sea quark 
mass: see figure^ The values of these curves at o" = are also given in table 101 

The extrapolation of the full measured ratios, including 0{a) contributions, 
requires an extrapolation of rjs as a function of a which has, of course, a considerable 
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Table 6: Results for the parameters of the sea quark mass depen- 
dence. Quantities directly used in the fitting procedure are in bold 
face. 



Lr45 


4.34(28) ■ 10-3 


Rf{a = 0) 


0.415(19) 


Wr45 


1.1(1.4) ■ 10-3 








-9.1(6.4) ■ 10-5 


Rn{a = 0) 


1.025(17) 


Wi?6845 


-5.52(48) ■ 10-3 







uncertainty. The behavior of the extrapolated curve is especially sensitive to the 
assumed form of the r/5-extrapolation for Rnss near zero. For instance, if the 
magnitude of the 0{a) contribution given by ps = ijsxs is finite at zero, which is 
reasonable to assume, then Rnss = ^ singularity near zero. 

This is a manifestation of the fact that different definitions of the "critical line" in 
the (/?, K)-plane, for instance by = or m^^^^ = 0, in general differ by lattice 
artifacts (in our case by 0{a)). If, however, r/s = ps/xs would have a finite value 
at cr = then there would be no such singularity. The two extrapolations shown in 
the lower part of figure |1] are examples of these two cases. 

Concerning the results on the parameters obtained from the sea quark mass 
dependence (table IHl and the second half of table one has to remark that the 
assumption of a quark mass independent lattice spacing a has an important effect 
on them. Assuming a quark mass independent Sommer scale parameter rg would 
change these results substantially. (There would be small changes in the first half 
of table IHl due to the somewhat different values of the quark mass ratios ai^2, too.) 
For instance, the values of A4//0 and A3//0 would come out to be 16.1(1.1) and 
30.4(2.9), respectively, instead of the values given in the tables. As it has been 
discussed above, the choice of a quark mass independent renormalization scheme 
requires a quark mass independent lattice spacing and is not consistent with a quark 
mass independent rg [251 128j . Nevertheless, it is plausible that in the continuum 
limit and in the limit of very small sea quark masses rg/a becomes independent 
from the sea quark mass and the differences between the values for constant tq and 
a disappear. 
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5 Summary and discussion 



The results obtained in this paper for the Gasser-Leutwyler constants (see tablesHHSl 
and|ni) can only be taken as estimates of the values in continuum. In order to deduce 
continuum values with controlled error estimates the left out lattice artifacts have 
to be removed by performing simulations at increasing /3 values and extrapolating 
the results to a = 0. Reasonable next steps would be to tune the lattice spacing 
to a ~ 0.13 fm on 24^ • 48 and a ~ 0.10 fm on 32'^ • 64 lattices. This would require 
with the TSMB algorithm by a factor of about 10 and 100 more computer time, 
respectively. Our calculations near a ~ 0.20 fm should be improved by going from 
16^ to 16^ • 32 lattices in order to improve the extraction of the physical quantities 
of interest. The number of sea quark masses considered should be increased to 5-6 
towards smaller values. This will decrease the overall statistical errors considerably. 
We hope to reach sea quark masses about m^ea — \rns on 16'^ • 32 lattices in the 
near future. 

General conclusions of the present work are: 

• Compensating 0{a) effects in the pseudo-Goldstone boson sector by introduc- 
ing 0{a) terms in the PQCh-Lagrangian itself is a viable alternative to the 
C'(a)-improvement of the lattice action. An extension to also treat 0{a?) ef- 
fects in the PQCh-Lagrangian is possible |2()| l21j and has been partially taken 
into account also in the present paper. 

• The observed 0{a) contributions in the pseudo-Goldstone boson sector are 
surprisingly small. The ratios of the 0{a) parameters in the NLO PQCh- 
Lagrangian to the quark masses r]s = Ps/xs ^-re in our present range of quark 
masses {^rus < nigea ^ f^T-s) at the few percent level. 

• Taking ratios of pion mass squares and pion decay constants at fixed gauge 
coupling and varying quark masses has the advantage that the Z-factors of 
multiplicative renormalization as well as all sorts of quark mass independent 
corrections cancel. 

• NNLO contributions in PQChPT are in our present sea quark mass range 
rather important. In fact, they are more important than the 0{a) lattice 
artifacts. This introduces new parameters in the multi-parameter fits which 
makes the fitting procedure more difficult. The situation will be better at 
smaller sea quark masses where the importance of NNLO terms diminishes. 

The present results strengthen the observation already made in our previous 
paper [Hj that the expected behavior dictated by PQChPT sets in rather early - 
at relatively large lattice spacings - once the quark masses are small enough. Our 
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present cut-off ~ 1 GeV is already a "high energy scale" from the point of view 
of the pion dynamics. As a consequence, it seems to us that the numerical study of 
the pseudo-Goldstone boson sector of QCD is perhaps the easiest field for obtaining 
new quantitative results about hadron physics by lattice simulations. 
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Figure 1: {RRn — 1), [Rnvv — ^) and {Rnys — ^) for the three differ- 
ent sea quark mass values (full lines). Beside the fit the unphysical 
contribution (proportional to rjs) is separately shown (broken lines). 
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Valence quark mass dependence of RRf 

le'* lattice Xh = 33.5(2.4) 
|i=5.10 = 1.0 I 

K=0.1760 T| = 0.07 



0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 
Valence quark mass dependence of RfW 



le'* lattice ):r = 33.5(2.4) 
p=5.10 = 1.0 

ic=0.1760 Ti = 0.07 




0.6 0.8 1 1.2 1.4 1.6 1.8 2 
Valence quark mass dependence of RfVS 



16* lattice 


Xr = 33.5(2.4) 




0=1.0 


ic=0.1760 


Tl = 0.07 








fit 



Valence quark mass dependence of RRf 
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Figure 2: {RRf-1), {Rfw-l) and {Rfvs-l) for the three differ- 
ent sea quark mass values (full lines). Beside the fit the unphysical 
contribution (proportional to rjs) is separately shown (broken lines). 
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Valence quark mass dependence of RRn 
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Figure 3: NNLO tree-graph contribution at k,2 = 0.1770 where the 
sea quark mass is given by Mr ~ 1 (broken lines). The full lines 
represent the total fits shown also in figuresUIS^ which are the sums 
of the continuum NLO, the 0{a) and NNLO terms. 



21 



Sea quark mass dependence of Rf 
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Figure 4: Rfss (^^'d Rnss- The full lines show the estimate of 
continuum contributions - without 0{a) terms. The broken lines 
near zero show the unphysical contributions (proportional to rjs). 
The other lines are different extrapolations of the measured values 
including 0{a) terms (see text). 
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